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Q : Abstract 



Oh. It is shown using Monte Carlo simulation that for low multiplicity 

events the single-event factorial moments are saturated by the sta- 



tistical fluctuations. The diverse of the event-space moments C Pt 



of single-event moments with the diminishing of phase space scale, 
called "erraticity" , observed in experiment can readily be reproduced 
by a flat probability distribution with only statistical fluctuations and 
therefore does not indicate the existence of chaos as suggested. The 
possibility of studying chaos in high multiplicity events using erraticity 
analysis is discussed. 
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Since the experimental observation in 1983 of unexpectedly large local fluctuations in a 
single event of very high multiplicity recorded by the JACEE collaboration [|IJ, the investiga- 
tion of non-linear phenomena in high energy collisions has extracted much attention || . The 
anomalous scaling of factorial moments [[J averaged over event sample called intermittency 
(or fractal) has been proposed for this purpose. Such kind of anomalous scaling has recently 
been observed successfully in various experiments [[]] || . 

Beside the moments averaged over event sample, the importance of the fluctuation of 
single-event moments inside an event sample has also been stressed ||. It is shown that this 
kind of fluctuation is related to the chaotic behavior of the system ||. A quantity /i 9 called 
entropy index has been introduced |7] as an adequate parameter in measuring the chaotic 
behavior. The positivity of entropy index fi q > is proved to be a criterion for chaos [Q. 
This method has been given the name of "erraticity analysis" H]. 

The idea of studying chaos through the fluctuation of single-event moments in the event 
space, in addition to the conventional study of intermittency (fractality) through the mo- 
ments averaged over event sample, is enlightening. However, the method of eliminating 
statistical fluctuations using factorial moments, which worked well for the moments aver- 
aged over event sample ||, can not be simply extended to the case of single-event moments, 
cf. Appendix. 

In this letter we will take this problem into account. We will show, using Monte Carlo 
simulation, that the single-event factorial moments are saturated by the statistical fluctua- 
tions when the multiplicity is low. The phenomena observed in experiments ]ID) , which were 
interpreted as a signal of chaos can readily be reproduced by a flat probability distribution 
with only statistical fluctuations. The idea || of studying chaos in high energy collisions 
using the distribution width of single-event moments in event space is meaningful only for 
high multiplicity events. 



In order to study the influence of statistical fluctuations on the chaos-study using single- 
event moments, let us rephrase the formalism || in terms of both the factorial moments and 
the probability moments parallelly. 

The factorial moment and probability moment Of™ for each event are defined as 

M 

jf E n m (n m - 1) • • • (n m -q + 1) 

\ m=l J 

M 

C (e) = M ^J2(pt ) y (2) 
i=l 

respectively. 

The moments F^ e ' (or C<p) may fluctuate greatly from event to event. In a sample 
consisting of a large number J\f of events, we get a distribution P \F^J (or P^C^^, 
which is normalized to unity. By taking the normalized moments of these distributions in 
event-space defined as 

C£ = (F^) /{F^Y , C$ = (C^ P ) /((%>)> , (3) 

we have a quantification of the fluctuation of the spatial patterns, i.e. we can investigate 
the full shape of the distribution and, especially, the way it changes with the resolution 
5 = A/M. The value of p in eqn.(3) can take any positive real number. If C^J(M) (or 
C^(M)) has a power law behaviour in M, i.e. 

C<5 (M) oc M^ F) (p) , CjJ (M) oc M^ C) (p) , (4) 

then the entropy index can be defined as, (For simplicity in notation, the superscript (F) 
and (C) will be omitted in the following.) 

(5) 



P =i 

It is easy to see that finite, nonvanishing positive values of fi q corresponds to wide dis- 
tribution of F^ (or C^), which in turn means unpredictable spatial pattern from event to 



event. By applying the measure to known classical chaotic system, it has been shown 
that fi q can be used as a measure of chaos in problems where only the spatial patterns can 
be observed and the positivity of fi q is a criterion for chaos. 
An alternative way of calculating \i q || is to express C PtQ as 

C M =(^ e)P >, (6) 



in which, 



With the definition 



*« = CjO/(Cf). (7) 



E fl = <<&«ln<&«), (8) 



we can obtain 



in the scaling region, i.e. where E g exhibits a linear dependence on InM. 

In order to study the influence of statistical fluctuations let us turn now to the application 
of this formalism to a simple model of dynamical fluctuations — the random cascading a- 
model and consider both the case of pure probability distribution without any particle 



distributed in the space and the case with particle distribution. 

In the random cascading a-model, the M divisions of a phase space region A are made 
in steps. At the first step, it is divided into two equal parts; at the second step, each part in 
the first step is further divided into two equal parts, and so on. The steps are repeated until 
M = AY/Sy = 2 U . How particles are distributed from step-to-step between the two parts 
of a given phase space cell is defined by independent random variable u) V j v , where j v is the 
position of the sub-cell (1 < j v < 2 U ) and v is the number of steps. It is given by fll |: 

u v ,2j-i = ^(1 + qt) ; L0 u>2j = i(l - or), j = l,...,2 u ~ 1 (10) 



where, r is a random number distributed uniformly in the interval [—1, 1]. a is a positive 
number less than unity, which determines the region of the random variable to and describes 
the strength of dynamical fluctuations in the model. After v steps, the probability in the 
mth window (m = 1, . . . , M) is p m = u\j 1 U2j 2 ■ ■ ■ Uvj u - A certain number N of particles will 
then be put in the M sub-cells according to the multinomial (Bernouli) distribution 

AH ^ 

The factorial moment and probability moment in each event can now be calculated 
according to eqn.(l) and (2) respectively. Then the moment C PA and entropy index \x q of 
the sample can be obtaioned using eqn.(3) and eqn.(9). 

Firstly, let us consider the case of pure probability distribution without any particle. It 



can be shown [|T^] that when the model parameter a is fixed to a constant value there will 
not be any chaotic behaviour, i.e. the entropy index fi will vanish. A positive fi together 
with the anomalous scaling of averaged moment (intermittency) can be obtained when a is 
distributed over a certain range. For simplicity let a be a random variable having a Gaussian 
distribution. The mean and variance of the Gaussian are both chosen as 0.22. The results 
of lnC*2, lnC Pi 2 and £ 2 as function of InM are shown in Fig. 1(a). It can be seen from the 
figure that both the intermittency (power-law behavior of the averaged moment C 2 ) and the 
chaos (diverse of C p ^ for the increasing of M) are reproduced. The intermittency index 2 
and entropy index /i 2 obtained through linear fit of lnC*2 and S 2 vs. InM (in the latter case 
only the last 3 points are used in the fit) are: 

02 = 0.056 , H2 = 0.016 (12) 

respectively. 

Now let us put in a certain number, say 9, of particles according to the Bernouli dis- 
tribution, cf. eqn.(ll). The results of lnF 2 , lnC P 2 and S 2 as function of InM are shown 



in Fig. 1(6). It can be seen from the figure that both the intermittency and the chaos are 
survived. However, the scale of C p ^ increases for an order of magnitude, showing a much 
stronger diverse of C Pi 2 with the increasing of M. The intermittency index <p2 and entropy 
index \i2 obtained in the same way as do in the previous no-particle case are 

2 = 0.058 , /i 2 = 0.41 (13) 

respectively. 

Comparing eqn.(13) with eqn.(12) we see that the anomalous scaling of the averaged 
factorial moment F2 is very similar to that of the averaged probability moment Ci- The 
corresponding intermittency indices are almost equal. This means that in the moments av- 
eraged over event sample the statistical fluctuations are very well eliminated by the factorial 
moment method pi]. On the contrary, the entropy index [1% for the 9-particle case is much 
bigger than that for the pure probability case, bigger for order of magnitude. This large 
"entropy index" is clearly comes from statistical fluctuations. 

In order to convince us of the dominant role of statistical fluctuations let us consider the 
case with perfectly flat probabilty distribution, i.e. p m = 1/M (m = 1, . . . , M), without any 
dynamical fluctuation. In this case, calculating with probability moments, both the averaged 
C q and the C p>q are identically equal to unity, and the lnC 9 , lnC Pj9 and S g all vanish, i.e. no 
intermittency, no chaos. 

After putting in 9 particles, the averaged factorial moments keep constant approximately, 
independent on M, cf. the upper figure of Fig. 1(c). (Its value moving down from zero is due 
to the particular normalization used.) Therefore, the intermittency index almost vanishes 
as expected. However, the C v ,2 diverse strongly, cf. the middle figure of Fig. 1(c). The 
intermittency and entropy indices obtained through linear fit are 

2 = 0.00046 , & = 0.49 (14) 



respectively. This large value of /i2 is totally due to statistical fluctuations. Note that 
in getting this result the random cascading a model has not been used and therefore the 
conclusion is model-independent. 

The results of flat probability distribution with 9 particles are compared with those from 



the NA27 data [10J in Fig. 2. It can be seen from the figure that the fit is very well, which 
means that the phenomena observed in this experiment, being dominated by statistical 
fluctuations, does not indicate the existence of chaos. 

The particle number 9 used above is about the average multiplicity of the moderate- 
energy hadron-hadron collision experiments, such as NA27 and NA22. When the number of 
particle increases, the influence of statistical fluctuations decreases. In Fig. 3 are shown the 
"entropy indices" fi q of flat probability distribution with different number N of particles. 
It can be seen that fi q decreases with the increasing of N and tends to zero when N goes 
larger and larger. At the upper-right corner of Fig. 3 are plotted the fits of fi q vs. iV to the 
empirical formula 

H q = A q e- b * N0 \ (15) 

The fitting parameters are listed in Table I. 

In Fig. 4 are shown the dependence of second order entropy index /i 2 on the multiplicity 
for both the flat probability distribution and the a model with Gaussian distributed a. It 
can be seen from the figure that when multiplicity is low, e.g. lower than 50, the two cases 
are almost undistinguishable. It is worthwhile noticing that in this case the 'entropy index' 
\i-2 for flat probability distribution is even higher than that for Gaussian distributed a. This 
means that the existence of chaotic behavior cannot be revealed using this method. When 
multiplicty is high, e.g. higher than 200, the entropy index fi 2 for Gaussian distributed a 
becomes greater than that for flat probability distribution and approaches the physically 
meaningful entropy index of probability distribution (solid line of Fig.4). 



It is interesting to try this kind of analysis with the central collision data of relativistic 
heavy ion collision experiments, where the multiplicities are very high. If the entropy indices 
fi q obtained in these experiments turn out to be considerably bigger than the value for flat 
probability distribution, shown in Fig. 3 and the full circles of Fig. 4, then that will be a signal 
of some interesting new physics. 

Table I Fitting parameters of /i g vs. N to eqn.(15) 



q 


\ 


\ 


2 


2.124 


3.954 


3 


2.397 


3.463 


4 


2.526 


3.072 



In summary, it is shown in this letter using Monte Carlo simulation that the single-event 
factorial moments are saturated by the statistical fluctuations. The diverse of the event-space 
moments C PA of single-event moments with the diminishing of phase space scale (erraticity) 
observed in experiment can readily be reproduced by a flat probability distribution with 
only statistical fluctuations and therefore does indicate the existence of chaos as suggested. 
Using "erraticity analysis" to study chaos in high energy collisions is meaningful only for 
high multiplicity events. 
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Appendix On the elimination of statistical fluctuations in the vertically and 
horizontally averaged factorial moments 



Firstly, let us recall the elimination of statistical fluctuations in the vertically averaged 
FM briefly. The basic assumption was that the statistical fluctuations are Poissonian: 

Qn = / -e~ s D(s)ds, (Al) 
Jo n\ 

where Q n is the multiplicity distribution in the ith bin, D(s) the dynamical distribution in 
the same bin. In eqn.(Al) we have used the notation n with no subscript to denote the 
number of particles in the zth bin || The multiplicity distribution Q n and the dynamical 
distribution D(s) are restricted by the normalization conditions 

oo 

n=0 
p CO 

/ D(s)ds = 1. (A3) 
Jo 

In order to check whether the equation (Al) is consistent with the conditions (A2) and (A3), 
insert (Al) into (A2) 

n=0 J[J \n=0 "" J 

Using the normalization of Poisson distribution 

00 n 

E ^ e ~ s = 1, W 

we get 



n=o n - 



X POO 

E Qn = / D( S )ds 



consistent with eqn.(A2) and (A3). 

It should be noticed, however, that the normalization (A5) of Poisson distribution is 
valid when and only when the summation for n is over all the positive integer (and 0) values 



of n. If some terms in the summation are missing, or if the summation is ended at some 
value n max of n, then it will be no longer equal to unity. 

Let us consider the elimination of statistical fluctuations in the 2nd order FM as exam- 
ple. Consider only the numerator f 2 = = n{n — 1) in the definition eqn.(l), and neglect 
the denorminator, which is only for normalization. We have 

yV e=l 

The crucial point is that, when the number of event Af in a sample is big, the bin- multiplicity 
nf 1 (= n) for fixed i can take value from zero to very large (constrainted only by energy 
conservation). Therefore, we have 

oo oo n r00 oo n 

</2>v = E n(n-l)Q„ = E / n(n-l)-e- s D(s)ds = / E TT" ^ S D(s)ds . (A7) 

n=0 n=0 J0 U - J0 n=0 \ U Z )' 

Let n — 2 = n' and make use of the normalization (A5) of Poisson distribution we get 

roo I 00 q n ' \ roo 

h = (/ 2 > v = / E V' S s2D ( s ) ds = / s 2 D(s)ds = 02. (A8) 

Jo \ n '=0 n ■ / Jo 

This is the well known result - - (vertically averaged) factorial moment (FM) equal to 
dynamical moment (DM), i.e. the statistical fluctuations have been eliminated successfully 
in vertically averaging over event sample. 

Let us check whether this conclusion is applicable also to horizontally averaging. There 
are 2 points that are different to the previous vertically averaging case. 

1) The number M of bins, unlike the number J\f of events, cannot be very large, may be 
only 1, 2, . . . , or at most some tens or hudreds. 

2) Even more important, the total multiplicity in a single event is always fixed. 
We can of course define a bin-multiplicity distribution of a single event as the 

distribution of bin- multiplicity n(= nf^) for fixed e in the M bins, but Qffi is nonzero only 
for at most M' different values of n, where M' < M. 

^ only for n — n x , n 2 , . . . , n M >- 



Among the M' integers ri\, 112, ■ ■ ■ , Um', there is certainly a maximum one n max . 

= when n > n max . (A9) 
Therefore, for the single event bin-multiplicity distribution, 



00 roo /"-max „n \ 

n=0 J ° \n=0 U - / 



Since the summation over n in the round bracket ends at a finite number n max and has only 
M' terms, it does not equal to unity and 

00 />oo 

/ D(s)ds. (A10) 
„=o 7o 

This means that eqn.(Al) is inconsistent with eqn.(A2) and (A3) for Q$. In other words, 
for the single-event bin-multiplicity distribution Qffl normalized as eqn.(A2) no normalizable 
D(s) which satisfies eqn.(Al) can exist. Therefore, the proof for the elimination of statistical 
fluctuations is no longer valid in the horizontal averaging case. 
State alternatively, if we write 

00 

(/ 2 > h = E^ e) ^-i), (^12) 

71=1 

then due to eqn.(A9) the summation in n doesnot really extend to infinity. If we neglect the 
non-existence of D(s) for a while and try to do the calculation similar to eqn.(A7),(A8), we 
will get 

(/ 2 >h = E^"l)Qi e) =E / n(n-l)-e- s D(s)ds 

n=l n=l J n] 

pOC f n ma x &n \ «oO / n max — 2 n' \ 



The normalization condition (A5) cannot be used here, and therefore 

(/ 2 )h ^ / s 2 D(s)ds. 
Jo 



(A13) 
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Figure Captions 



Fig.l Logarithm of averaged moment I11C2 (lni^), lnC Pj 2 and £2 as function of 
InM for (a) the probability moments of a model with Gaussian-distributed 
a; (b) the factorial moments of the same model with 9 particles put in; (c) 
the factorial moments of a fat probability distribution with 9 particles put 
in. solid lines are linear fit. Dashed lines are for guiding the eye. 

Fig. 2 InCp^ (q = 2-4) vs. InM for a flat distribution with 9 particles put in 
(full squares) compared with the experimental results of NA27 data (open 
circles). Data taken from Ref.fllU . 



Fig. 3 The dependence of the entropy indices fi q on the number of particles. The 
solid lines are the results from the empirical formula (15) and Table I. 

Fig. 4 The dependence on the number of particles of the entropy indices H2 calcu- 
lated from factorial moments. Full circles - - flat probability distribution 
with only statistical fluctuations; open triangles - - Gaussian distributed 
a. The solid line is the result from probability moments with Gaussian 
distributed a. The dashed lines are for guiding the eye. 




Fig.l Logarithm of averaged moment I11C2 (lnF 2 ) and lnC Pi 2 as function of InM 
for (a) the probability moments of a model with Gaussian-distributed a; 
(6) the factorial moments of the same model with 9 particles put in; (c) the 
factorial moments of a fat probability distribution with 9 particles put in. 
Solid lines are linear fit. Dashed lines are for guiding the eye. 
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Fig. 2 lnC Pt g (q = 2-4) vs. InM for a flat distribution with 9 particles put in 
(full squares) compared with the experimental results of NA27 data 
(open circles). Data taken from Ref . |T0| . 




Fig. 3 The dependence of the entropy indices jj q on the number of particles. 
The solid lines are the results from the empirical formula (15) and 
Table I. 
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Fig. 4 The dependence on the number of particles of the entropy indices ^ 
calculated from factorial moments. Full circles - - flat probabil- 
ity distribution with only statistical fluctuations; open triangles - 
Gaussian distributed a. The solid line is the result from probabil- 
ity moments with Gaussian distributed a. The dashed lines are for 
guiding the eye. 



